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Abstract 



In this short note, we describe the preparation of updated templates for the interpretation 

of SUSY results from the LHC in the context of mSUGRA. The standard {mo,m,i/2) plane is 
shown for fixed fi > and nit = 173.2 GeV. Two scenarios are considered: (1) Aq = GeV 
and tan(/3) = 10 and (2) Aq = —500 GeV and tan(/3) = 40. In each case, the universal scalar 
mass parameter mo varies in the range mo G [40, 3000] GeV, while the universal gaugino mass 
parameter mi/2 varies in the range mi/2 € [100, 1000] GeV. We delineate notable regions in 
parameter space, including the region with a charged LSP (stau), the LEP2 reach, and the 
cosmologically preferred region with 100% neutralino dark matter. The templates also show 
mass contours for a few key particles (gluino, squark and Higgs boson). The mass spectrum is 
calculated with the SoftSusy-3.2.4 package, while the neutralino relic density is obtained with 
MicrOMEGAs version 2.4. 



I. INTRODUCTION 

Low energy supersymmetry (SUSY) [IHTj is a primary target of the LHC collaborations 
in their quest for new physics at the TeV scale. Unfortunately, even the most minimal 
supersymmetric extension of the Standard Model comes with a large number of a priori 
unknown parameters (superpartner masses, mixing angles and CP violating phases). The 
experimental exploration of such a large parameter space is impractical. Therefore, it 
has become customary to present experimental limits from searches for supersymmetry in 
terms of simple benchmark models with very few input parameters, see e.g. [SHT^. The 
most popular SUSY benchmark model is the minimal supergravity (mSUGRA) model [13]- 
[15], a.k.a. the "constrained MSSM" (cMSSM). In this model the SUSY mass parameters 
are input at the grand unification (GUT) scale, conventionally defined as the scale where 
the two gauge couplings gi and g2 meet. At the GUT scale, all scalars in the model 
have a common mass mo, while all gaugino fermion superpartners have a common mass 
mi/2- The model has two additional continuous parameters: Aq, a common trilinear scalar 
coupling at the GUT scale, and tan(/3), the ratio of the two Higgs vacuum expectation 
values, plus a discrete parameter, the sign of the higgsino mass parameter /i. In order to 
obtain the physical SUSY mass spectrum, the parameters mg, mi/2 and Aq are evolved via 
the renormalization group equations (RGEs) from the GUT scale down to the electroweak 
scale, where the radiatively corrected mass spectrum is computed. Since the gauge and 
Yukawa couplings of the Standard Model are input at the electroweak scale, the procedure 
must be iterated until it converges on a stable solution. There are several state of the art, 
publicly available codes on the market which can do these calculations, and they generally 
give similar results in the bulk of the parameter space [16]. Here we use the SoftSusy-3.2.4 
software package [T7] . 

TABLE 1. Relevant input parameters for SoftSusy mass spectrum calculations 



Parameter 


Value 


Oisimz) 


0.1184 




127.934 


nit 


173.2 GeV 


nib 


4.2 GeV 


sign(;u) 




+ 


Aq 


GeV 


-500 GeV 


tan(/3) 


10 


40 



In this note we provide details of the design of mSUGRA templates which can be 
used for interpretation of supersymmetry searches at the LHC. Given the ubiquitousness 
of the mSUGRA model, it has become customary to present exclusion limits from SUSY 
searches in the mo—mi/2 parameter plane, for fixed values of the remaining mSUGRA (and 
Standard Model) parameters (shown in Tablejl]). In anticipation of the improved reach with 
2011 and 2012 data, we extend the previously considered ranges of mo and mi/2, allowing 
them to vary in mo G [40,3000] GeV and mi/2 G [100, 1000] GeV, correspondingly. Two 
different mo — mi/2 slices are considered. The first has a moderate value of tan(/3) = 10, 
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iTio [GeV] trio [GeV] 

(a) tan(;9) = 10 (vl) (b) tan(;9) = 10 (v2) 



FIG. 1. Template versions 1 (left) and 2 (right) for the tan(/3) = 10 scan. The features are 
described in the text. 
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FIG. 2. Template versions 1 (left) and 2 (right) for the tan(/3) = 40 scan. The features are 
described in the text. 



while Aq is taken to be = 0. The second mo — mi/2 shce has a higher value of 
tan(/3) = 40 and a different Aq = —500 GeV. In each case, the sign of the fi parameter 
is positive, and the top mass is taken to be rrit = 173.2 GeV. In what follows, we refer to 
these two mo — mi/2 scans by their respective tan(/3) parameters. 



II. MSUGRA TEMPLATES 



We begin by presenting the templates in their final form before discussing the details 
of their respective features. The templates are displayed in Figs. [T] and [2] for tan(/3) = 10 
and tan(/3) = 40, respectively. The calculations are performed across an niQ — mi/2 grid. 
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The LHC mSUGRA scans [181 sample this grid in 20 GeV intervals for both mo and 
mi/2 with about 10k events for each point. A finer granularity is of course desirable 
but is computationally prohibitive, as each grid point has to be propagated through the 
generation, simulation, and reconstruction chain. The storage space required for hosting 
these datasets would also grow quadratically with finer granularity. Since our goal is to 
simply provide a canvas which delineates theoretical and experimental regions of interest 
in the mg — mi/2 plane, we only need to make use of the mass spectrum calculations, 
and these come in concise SLHA files from the SoftSusy output. This allows us to pursue 
a much finer granularity of 5 GeV intervals. In the high tan(/3) scan, grid points with 
mo < 220 GeV do not yield a suitable neutral dark matter candidate, thus the scan 
effectively begins at mo = 220 GeV. 



III. CONSTRAINTS ON THE MSUGRA PARAMETER SPACE 

Within the mSUGRA parameter space considered here, there are regions which are 
not considered to be viable for hypothesis testing, since they are subject to constraints 
imposed by theory or experiment. 



A. Theoretical constraints 



1. Charged LSP 



One nice feature of the mSUGRA model is that throughout most of the mg — m.1/2 pa- 
rameter space, the lightest supersymmetric particle (LSP) is a neutralino, which represents 
a potential dark matter candidate. If a suitable discrete symmetry (such as i?-parity) is 
imposed on the model, the neutralino does become a dark matter candidate and its relic 
density can be readily computed from the thermal freeze-out (see Sec. Ill B 3 below). More 
importantly for collider searches though is the fact that neutralinos are very weakly in- 
teracting particles, and do not interact inside the detector. Thus mSUGRA parameter 
space points with neutralino LSP will be relevant for SUSY searches utilizing the missing 
transverse energy in the event. 

As can be seen from Figs. [T] and |2| though, at low values of mo there is a small corner 
in the rriQ — mi/2 plane (colored dark gray and denoted "f = LSP"), where the lightest 
particle in the SUSY spectrum is a r slepton (stau f). Being an electrically charged 
particle, the stau is not a good dark matter candidate, therefore the dark gray regions in 
Figs. [1] and [2] are cosmologically disfavorecQ 



2. Correct electroweak symmetry breaking 

Another important theoretical requirement on the model is that the global minimum of 
the scalar potential should break the electroweak symmetry, leaving QED and QCD intact. 

^ This conclusion can be avoided if a) the stau is not the true LSP in the model, but slowly decays 

to a superWIMP such as the gravitino [12] or b) if i?-parity is violated and the stau decays to SM 

particles, with no sizeable missing transverse energy in the event. In either case, the corresponding 

collider signatures fall into the scope of the exotica group. 
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Within the mQ — mi/2 plane one does find regions with the wrong gauge symmetry breaking 
pattern. For example, at very large values of mo, one finds an area with no electroweak 
symmetry breaking, as indicated by an unphysical value of /i^ < 0. In Figs, [l] and [2] 
that region is colored light gray and denoted "No EWSB" . (At very small values of mo at 
large tan(/3), one may encounter an area with a charge-breaking vacuum, as signalled by a 
negative stau mass squared (M| < 0) [20]. For simplicity, we do not specifically delineate 
that region, since it is already well inside the stau LSP region.) 



3. Convergence on reliable solutions to the RGEs 

As already mentioned in the Introduction, the numerical programs used to solve the 
MSSM RGEs use an iterative procedure, since some inputs (the gauge and Yukawa cou- 
plings and tan(/3)) are specified at the weak scale, while others (mo, mi/2 and Aq) are 
defined at the unification scale. Given that the RGB's are very nonlinear functions, one 
might expect that in some situations the iterative procedure will not converge and instead 
will exhibit chaotic behavior. Indeed, this turns out to be the case in a narrow strip 
(colored light blue), which is adjacent to the "No EWSB" region in Figs. [l| and g 

Like all other SUSY spectrum calculators, SoftSusy checks for convergence after each 
iteration, and eventually terminates after a given fixed number of iteration^ One should 
keep in mind that even though SoftSusy warns the user of the non-convergence problem. 



it still provides a mass spectrum, albeit an unreliable one. This is illustrated in Figs. 3 (a 
and 3(b) I which show the masses of the two charginos obtained in SoftSusy as a function 



of mo for a horizontal slice across the mo — mi/2 planes of Figs. [T] and [2] at a fixed value 
of mi/2 = 220 GeV. Here, the region of RGE non-convergence as reported by SoftSusy is 
shaded in light blue. In this region, the expected value of as calculated from minimizing 
the radiatively corrected effective potential, is rather small, and the lightest chargino xt 
is higgsino-like. The oscillatory behavior of its mass M-+ at high mo is an indication of 
the convergence problem. Since the value of in the non-convergence region is small, 
this region is always found to buffer the gray "No EWSB" region in which /i^ < 0. It is 
important to emphasize that non-convergence of the RGEs by itself does not necessarily 
mean that such points are ruled out or disfavorec^^ — they may very well represent physi- 
cally viable models and one should not mistakenly incorporate them into the "No EWSB" 
region where /i^ becomes negative. Those are indeed two separate regions and deserve to 
be distinguished on the canvas. 



In both Figs. 3(a) and 3(b) we observe perfect agreement between the chargino masses 
as calculated by our private SoftSusy scan (solid lines, 5 GeV granularity) and those 
calculated by the LHC scan [TS] (red markers, 20 GeV granularity). While for this slice 
of mi/2 = 220 GeV, the LHC scans do not extend into the regions where non-convergent 
RGE's are reported, it does happen in rare cases (a few out of a thousand) — the LHC scan 
would sample a point that falls just over the boundary of the region of non-convergence. 
We have checked that the chargino masses obtained in those cases do coincide with the 

^ We have verified that increasing the maximum number of iterations does not alleviate the convergence 

problem, i.e. the chaotic behavior is intrinsically present in the system. 
^ Non-convergence of the RGEs means simply that — the mass spectrum calculation is unreliable. In 

those cases, the SoftSusy manual [17j recommends using an extrapolation from the neighboring regions 

where the program did converge. 




No EWSS 
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m„ [GeV/c^] 

(a) tan(/3) = 10 




FIG. 3. Masses of the two charginos Xi and X2 GeV) as a function of uiq for mi^2 — 220 
GeV for the tan(/3) = 10 (top) and tan(/3) = 40 (bottom) scans. The non-convergent region is 
shaded in light blue, while the "No EWSB" region is shaded in gray. The green shaded region is 
ruled out by chargino searches at LEP. 



result that one would obtain by extrapolating from the convergent region, so no ill-effects 
are expected from this occasional trespassing. 



It should be noted that while all of these theoretical constraints are ascertained with 
the SoftSusy spectrum calculator, they are generally understood and expected constraints 
and appear when using other spectrum calculators as well (Suspect [2T|, SPheno [22] , 
IsaJet [23], etc.). 
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B. Experimental Constraints 



Results from previous experiments can also be interpreted in the context of the 
mSUGRA parameter space shown in Figs. [T] and |2} 



1. Chargino mass limits from LEP2 




LEP2 searches for chargino pair production are probing the two areas of the mo — 
plane which exhibit a light chargino (at low mi/2, where Xi is wino-like, and at large 



mo where /i becomes small and Xi is higgsino-like). In Figs. [ljand[2jthe corresponding 

regions are colored green and can be understood as excluded at 95% C.L. by LEP2 ^23]. 

For reference, the lightest chargino mass M-+ is plotted across the mo — m.1/2 grid in 
I 1 I 1 •''-1 ' 

Figs. 4(a) and |4(b) for tan(/3) = 10 and tan(/3) = 40, respectively. 



2. Slepton mass limits from LEP2 

LEP2 has also performed searches for direct slepton production. Explicitly, the follow- 
ing exclusion limits are imposed [24j: 

m(gij) < 100 GeV A m(x?) < 85 GeV (1) 
m{fiR) < 95 GeV A m{j^R) - m(x?) > 5 GeV (2) 
m(fi) < 86 GeV A m(fi) - m(x?) > 7 GeV (3) 

The relatively small region which satisfies this condition can be found at low mo and low 
mi/2 and is colored yellow in Figs, [l] and [2j It can be understood as excluded at 95% C.L. 
by LEP2. 
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3. Dark matter relic density 



tan(B)=10, A = GeV, li > 0, m = 173.2 GeV tan(p)=40, A = -500 GeV, li > 0, m = 173.2 GeV 




mp [GeV] trio [GeV] 

(a) tan(/3) = 10 (b) tan(/3) = 40 

FIG. 5. Dark matter relic density VLh? as a function of rriQ and mi/2 for tan(/3) = 10 (top) and 
tan(/3) = 40 (bottom). 



In addition to collider constraints from direct searches for superpartners, one may also 
wish to consider the cosmological constraint on the dark matter relic density (which is 
indirect and perhaps more speculative, since it involves certain assumptions about the 
composition of the dark matter as well as the thermal history of the Universe). We calcu- 
late the neutralino relic density across the mo — "021/2 grid with MicrOMEGAs version 2.4 
[25] using SoftSusy-3.2.4 for mass spectrum calculations. The resulting O.-yh? is displayed 



in Figs. 5(a) and 5(b) for tan(/3) = 10 and tan(/3) = 40, respectively. One can see that 
throughout most of the parameter space, the neutralino relic density is quite larger than 
the WMAP preferred range [26] 



0.102 < VL^K^ < 0.123. (4) 

Model points satisfying eq. ^ yield the appropriate amount of dark matter to within 3a 
of the measured value of 0.1123 ± 0.0035 and are denoted by magenta filled circles in 
version 2 of the templates on Figs, and[2|^b). As expected, the proper rehc density is 
obtained in two areas: 1) the neutralino-stau coannihilation region [27] near the boundary 
of the "f LSP" region and 2) the "focus-point" region at large mo [281 EH], where the 
neutralino LSP picks up a non-negligible higgsino component [3Jlj. Other notable features 
on the plot in Fig. |5]^a) include the opening of the neutralino annihilation channel into 
pairs of top quarks (near mi/2 = 430 GeV) and pairs of Higgs and/or gauge bosons (in 
the neighborhood of mi/2 = 300 GeV). 



IV. MASS SPECTRUM 

It is a tradition to display a few meaningful isomass contours for particles of interest 
on the mSUGRA canvas. In the past, the gluino and squark masses were displayed in 
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250 GeV intervals]^ We continue to display these contours, opting for 500 GeV intervals 
to reduce congestion. In Version (2) of the templates we go further and provide two 
additional contours for a Higgs mass of 115 and 120 GeV, respectively. These values are 
particularly relevant at this stage as the results of the Higgs searches at the LHC may 
soon make contact with this parameter space [3Tj. As a reference. Figs. [6| [7| and [s] show 
the masses of the gluino, squark (average), and light CP-even Higgs boson, respectively, 
across the mo — mi/2 plane. 



tan(B)=10, A = GeV, li > 0, m = 173.2 GeV 

I ^ — Q * 
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tan(p)=40, = -500 GeV n > 0, = 173.2 GeV 




500 1000 



1500 2000 2500 3000 

m, [GeV] 



(b) tan(/3) = 40 



FIG. 6. Gluino mass in GeV as a function of mo and mi/2 for (a) tan(/3) = 10 and (b) tan(/3) = 40 



tan(B)=10, A„ = GeV u > 0, m = 173.2 GeV 



tan(p)=40, A I, = -500 GeV n > 0, = 173.2 GeV 
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(b) tan(/3) = 40 



FIG. 7. Average first generation squark mass in GeV as a function of niQ and for (a) 

tan(/3) = 10 and (b) tan(/3) = 40. 



Here we take the average of the four first generation squarks. These are observed to always be within 
8% of one another. 
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(b) tan(/3) = 40 

FIG. 8. Higgs mass in GeV as a function of uiq and mi/2 for (a) tan(/3) = 10 and (b) tan(/3) = 40. 



V. APPLICATION 

It is expected that Version (1) of the templates is used for public presentation. Both 
versions, however, are available as TCanvas objects and C-macros in the ROOT framework. 
These can be downloaded from the web at 

• tan(/3) = 10 

tier2.ihepa.ufl.edu/~remington/SUSY/mSUGRA/GridTanblO_vl.root (.C) 
tier2.ihepa.ufl.edu/~remington/SUSY/mSUGRA/GridTanblO_v2.root (.C) 

• tan(/3) = 40 

tier2.ihepa.ufl.edu/~remington/SUSY/mSUGRA/GridTanb40_vl.root (.C) 
tier2.ihepa.ufl.edu/~remington/SUSY/mSUGRA/GridTanb40_v2.root (.C) 

and are also available as ancillary files from the arXiv version of this note. 

The following code can be used in ROOT (v.27) to draw the canvas from the root files 
above: 

TFile f ("GridTanblO.vl.root' ; 

TCanvas *c = (TCanvas*) f .GetC 'GridCanvas") ; 

c->Draw() ; 

MyExclusionContour->Draw( ' 'SAME' ') ; 

As a final check, we survey the model points that are expected to be in the official 
LHC scans |lSj and we represent them symbolically with a "+" on top of our templates to 
assess the coverage. This can be seen in Figs, [ojand 10 for tan(/3) = 10 and tan(/3) = 40, 



respectively. For tan(/3) = 10 we see that there is full coverage everywhere. For tan(/3) 
40, however, we see a small envelope in between the dark matter preferred region and the 
non-convergent RGE region where the LHC scan does not sample. Our recommendation 
is that future scans should be extended to include this interesting region. Analysts should 
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be prepared to perform a small extrapolation in this region in the case that their exclusion 
contour terminates nearby. 
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<j 900 
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FIG. 9. LHC scan points (+) overlaying the mSUGRA template for tan(/3) = 10. 



^ 1000 i 
0) 

900 1 
S 8001 



iiiiiii i iiii i iiiiiiiiiliii i iii i iiiiiiiiiiiiiiliiii i iiiiiiiiiiiiii i iTrj;: 




500 1000 1500 2000 2500 3000 

mo [GeV] 



FIG. 10. LHC scan points (+) overlaying the mSUGRA template for tan(/3) = 40. 
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